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- ■ ABSTRACT 

Q ■ The ultraviolet asymptotic behavior of the gluon propagator is eval- 

. uated in the maximal Abelian gauge in the SU{2) gauge theory on the 

C . basis of the renormalization-group improved perturbation theory at the 

one-loop level. Square-root singularities obtained in the Euclidean do- 
lO ! main are attributed to artifacts of the one-loop approximation in the 

maximal Abelian gauge and the standard normalization condition for the 
CN [ propagator used in our study. It is argued that this gauge is essentially 

nonp erturbative . 
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1 Introduction 

Maximal Abelian Gauge (MAG) ilj offers one of the best instruments for theoretical studies 
of dynamics of the gauge fields in QCD at a large scale (~ 1 fm). In this gauge, the boson 



o • fields can be naturally divided into the "Abelian" and "non- Abelian" (or, in other words, 
"diagonal" and "off-diagonal") components. Assuming that the "diagonal" components give 
D ', a leading contribution to the path integral,^ the Yang-Mills theory on a lattice can be 
' reduced to the dual Abelian Higgs model. In so doing, confinement of quarks and gluons 
comes about from monopole condensation in the dual Abelian Higgs model. 
^ ■ Therefore, propagators of "diagonal" and "off-diagonal" gluons are of particular interest 
■ because, from a naive point of view, their ratio indicates whether the Abelian dominance 
hypothesis is true or not. Recently 13], the behavior of the "diagonal" and "off-diagonal" 
gluon propagators was estimated numerically in the framework of lattice gauge theory. At 
small momenta it was found that the propagator of the "Abelian" field components domi- 
nates. In this work, we study asymptotic behavior of the propagators of the " diagonal" and 
"off-diagonal" gluon at large momenta in the framework of perturbation theory in the con- 
tinuum field theory^. It should be noted that the interest in this gauge was quickened by the 
study ,5.^, where it was demonstrated that the SU{2) Yang-Mills theory can be renormalized 
in this gauge. 

-"^This assumption is referred to as the hypothesis of AbeHan dominance 

^We work in the Euclidean space, throughout this paper, = Po + Pi + P2 + pi, that is, Euchdean 
momenta are positive. 



1 



A ijrauge nxing 

In practical calculations in a lattice Yang-Mills theory, the maximal Abelian gauge is fixed 
by finding a stationary point of the functional 

$ = 5^Tr([/.,^a3t/V3) (1) 



with respect to the gauge transformations of the type 

A : f/,,^ At.f/,,^A,+^, (2) 

where Ux,/^ G SU{2) are the link variables and A^. is a SU (2)-valued function defined at sites 
of a lattice. In other words, the gauge condition is fixed by the relations 



5$ 

0. (3) 



SA 

At infinitesimal gauge transformations, 

A = exp{-zUa(ra/2) = cos — sm — , (4) 

the variation of the functional (0) takes the form 

(5$ = Tr [Ux^f,{a3 + Mi{x + il) + M2{x + jj))Ul^{cT3 + Mi{x) + M2{x))] 

where 

Miix) = -e3"^cu,(x)a,, M2ix) = i {5'''5'' - S'^'S^') UaU.a,. (5) 
Using the decomposition 

, -V^ f uo + ius U2 + iui\ 

''^ ^ V + iUl Uq - tU3 J 

n=l ^ ' 

we define the vector potentials by the formula -u" = —^A°^ (a = 1,2,3). In this case, 
the gauge condition has the form 

aWl{{ul-iul){ul-iul)) + 2ul{ul-iul) = 0, (6) 
«Vj((«J + mJ)(«i + my) + 2«J(«i + mJ) = 0, 



where Vf is the operator of backward derivative on a lattice: uUx) 

Having regard to the relation m° = 1 + 0{g^a?), we find that, in the leading order in the 
lattice spacing, the gauge condition ® has the form 

{d, + zgAl){Al + zAl) = Q. (7) 
{d,-igAl){Al-iAl) = Q. (8) 



Here and below, this naive limit of the above gauge conditions is named maximal Abelian 
gauge in the continuum theory. The continuum theory in this gauge can conveniently be 
quantized in the BRST formalism. To do this, we should introduce the gauge-fixing terms 
as follows: 

^^GF = -^ H + ^9Al)iAl + ^Al)\'-X {d,Al)\ (9) 

where a and /§ are the gauge parameters. The former term is responsible for the gauge 
condition (jZj), the latter term is needed to fix the remaining U{1) gauge arbitrariness. The 
case a = (3 = 1 corresponds to the case = 1 from article |4j. 

The BRST invariant Lagrangian of the Yang-Mills field in the maximal Abelian gauge 
has the form 
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+ ^ge^' [a^Ald^Al - a^A^d^A^ + A^A^d^a, - a^A^d^A^ 

+a^Ald,Al - AlA%a, - ^a^A^d^Al + ^a^A'^d^A^) 
- tge'''a,iC%C'' - C^'d^C') - tgKe'''C%iAlC'') 
+ ^g^S""' (^26^J,, - (1 - + V^-p)) AI^Ala^a^ 

+ \g\5,,5,„-5,,5,MlAlA';,Ai 

where the gluon field A"^ is divided into the "Abelian" and "non-Abelian" components as 
follows: 

where are the SU{2) generators; indices denoted by capital letters are assigned one of the 
integers 1,2,3 and indices denoted by small letters — 1, 2. The gluon fields with the indices 
denoted by small letters are named "off-diagonal" and the Abelian component defined 
by the equation (fTUj) is named the field of the "diagonal" gluon. Let us also introduce the 
notation: 

A detailed discussion of the parameter k can be found in here we set k = 1. 

3 Feynman rules 

The expressions for the gluon and ghost propagators and three- and four-particle vertices 
are readily obtained by standard techniques. The results are shown in the table 
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rropagators 


\a^[p)au[—p)) 


p2 \o^Jiu U PJ p2 j 




( S n n,\P^'P''\ 


{C\v)C\-v)) 


p2 


{C''{p)C\-p)) 


.gab 


Three-particle vertices 


{a^{p)A-{q)Ai{r)) 


-we'' ((? - r)^ + (r - p + f )/,^ 


{C\p)C\q)a,) 




{C'\p)C\q)Al) 


-nge^^p^. 


Four-particle vertices 




-gH-' [25^Jp„ - (1 - j;){5^p5,„ + 




„2 1 ^ab^cdj I ^ac^bdj _i_ ^ad^bcj \ 
g J-pv,p(T ~r t t J-pp^ua ~r t t J-ijl(t,i>p) 

IfjLU^pa i^^pp^ua ^pa^up) 




-i2gH'^'5^, 


{C^C'AlAt) 






g^ae^^e"'^ 



It should be noted that the "diagonal" ghost interacts with no field, whereas the 
"diagonal" anti-ghost interacts with some fields. Therefore, the "diagonal" ghosts should 
be disregarded in the loop expansion of perturbation theory. 

4 Gluon Propagator to One-Loop 

The diagrams contributing to the polarization operator of the "diagonal" gluon are shown 
in Fig. JH). 




Figure 1: Polarization operator of the "diagonal" gluon. 

Our computations are performed using the dimensional regularization techniques with 
the space-time dimension D = 4 — 2e; in so doing, the contribution of the tadpole diagrams 
vanishes. 

The contribution of the non-vanishing diagrams to the gluon polarization operator has 
the form 

n- = Tf^(Vp^-^.^^)(f^+(^ + 3« + Y))' (10) 
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where 



L = 7ij + In — — 

e \ 



and 'Je is the Euler constant. 

Performing analogous computations for the "off-diagonal" gluon (see diagrams in Fig. 
we arrive at 



Kip) 



+ 




Figure 2: Polarization operator of the "off-diagonal" gluon. 
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K(P) = - '^"'if^ (( W - P,P.)T"^^ip') + P,P.L""{P')) , (12) 

where 



..V) . ,-lI,,,|).,(-«4-^-|), (13) 



1 13 3 3\ l/^..o^3/5 5 
---F- ^U^H -7-h3/? 

2 a 2a a^/ 2a \ a a 

The expressions ()10|) and (fT^ involve both divergent and finite parts (the divergent part 
was computed in [H]). The finite part of the longitudinal component of the "off-diagonal" 
propagator p2|) involves terms proportional to ^ and which are singular when a — 0. 
However, they give no contribution to the expressions for physical quantities. Their presence 
in the expression for the propagator is due to the fact that the interaction vertex in this gauge 
is proportional to ^, where a is the gauge parameter that appears in the expression for the 
unperturbed propagator. For this reason, a — > corresponds to very dangerous limit and 
thus it is safe in the perturbation expansion to assume a nonzero value of a. 

In the MS subtraction scheme, the counterterms added to the Lagrangian used to com- 
pensate for the divergent terms are as follows: 

^C^ounter = ^^^^^ " ^ {d,Al - d^A^f + ^{O^Alf + (14) 



[d^Qy - d^a/,) + -^[o^a 
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where 



= 0, 



and the modified parameter e' is defined by the formula 

- = --7E + ln(47r). 

5 The renormalization group equation 

5.1 The case of Lorentz gauge 

First we consider the renormahzation group equation in the Lorentz gauge in order to demon- 
strate some features of its solution in the maximal Abelian gauge. The Lorentz gauge can 
be specified by the gauge fixing term in the Lagrangian as follows: 

A£gf = ^(^Kf- (16) 
The gluon propagator in this gauge is parametrized by only one scalar function 



(P) = ^ - ^) G-" {^^,, 9, a) + a Mil 1 . (17) 



The renormalization group equation has the form 



where 



/3(^) = /^^, 7(^) = ^^, %) = /^a;^- (19) 

The gauge parameter a is considered as yet another coupling constant |HI. 

In the dimensional regularization approach, the renormalization group functions f3{g) 
and 7 ((7) are determined from the relations 

^, , 1 / dai \ , , 1 dci . . 
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where ai and ci are the coefficients of - in the — expansion of the renormahzation factors 

e e 

Za and Zg-. 

Z, 



a 



i+-f:^ (21) 



^ n=0 



Za = l + V^. 

n=0 

In the one-loop approximation, the [3 function is independent of gauge and subtraction 
scheme. For the SU (2) gauge field it is given by 

In the Lorentz gauge, there are no contributions to the longitudinal component of the 
gluon polarization operator and there is no counterterms to the gauge fixing term in the La- 
grangian. From this observation we obtain the relation between the renormahzation factors, 
Z\ = Za, that determines the dependence of the gauge parameter a on the normalization 
point. In the case of the Landau gauge (a = 0) one can consider 6{g) = (see iQj); in the 
case a 7^ 0, the function 6{g) is not trivial. 

Here we solve the equation p8|l for 7 = cg'^,l3 = —bg^. First we note that, when 

1 

1(9) = 0, any function F(z) with z = — + 2b\n— provides a solution of the equation CHI)- 

9 

This function is specified by the initial (or boundary) conditions. That is, if the sought-for 
function is given over a line transversal to characteristic curves of the equation then the 

2 

function F{z) can be determined and thus the sought-for function is known over the (5', ^) 
plane. 

A useful method of solution of the equations of the type (fT^ with 'y{g) 7^ is based on 
a treatment of the solutions of the respective auxiliary equation 

where A = In — . Then the sought-for function G^°^(X, g, a) is determined from the condition 

l^(A,a,^?,G^n=0, (24) 

where V{X,a, g,G^'^^) is any solution of the auxiliary equation (|^. The solution of the 
auxiliary equation ()23|) is provided by any function which is constant along each characteristic 
curve. The characteristic curves are defined by the equation 

dX _ dg _ da _ dG'^'"' 

- ~ 2j{g,a)G-on- (25) 

Integration of these equations gives 

^ + = (26) 
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_ (IS) 13 - 3a 



3a 



a 



(27) 



(28) 



Each of these equations defines three-dimensional hypersurface and a characteristic curve is 
an intersection of these hypersurfaces. Thus the solution of the equation (f^H|) is given by 



LOR\ 



a 
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13 -3a 1 11 



3a ' c/2 247r2 
where v{x,y,z) is an arbitrary sufficiently smooth function. 



(29) 



50 

25 

13/3 



-25 



-50 



0, = 1 

f C2 = 0.01 

; C2 = -0.3 



Figure 3: Projection of characteristic curves onto the {a,g) plane. 



The condition ()24|) implies that 



^TOR ^/ ^(11)13 — 3a 1 11 
' ^ 3a ' ^2 ^ 247r2 



A 



(30) 



Formula ()30|) . where ^{y,z) is an arbitrary function, represents a general solution of the 
renormalization group equation in the Lorentz gauge. A specific form of the function $(?/, z) 
(a particular solution) is determined from the boundary conditions; it is natural to fix them 
at the hyperplane A = 0. The behavior of the characteristic curves makes it possible to 
determine the asymptotic behavior of the propagator for — >■ cx) and some fixed values of 
g and a provided that its behavior for g — >• 0, a — >• 13/3, and some fixed value of /i (say. 



/i^ = p^) is known. In the latter domain, perturbation theory works and thus the condition 
needed for the evaluation of the propagator at large momenta is fulfilled. In the leading 
order of perturbation theory, G = 1 and higher-order corrections are small. For A = 0, we 
obtain the equation 

(31) 



$(2/,z) = -, 
a 

which can be used for the determination of the dependence of $ on y and z. 

;^ 13-3a , 1 , , „ 13 



On this 



hyperplane, y = g 



3a 



and z = — and, therefore, 
9^ 



a 



3 1 



yz^ 



13 ^ 

22 y 
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Thus we arrive at 

Hy,z) = ^ [l + yzi-^^)y (32) 

Now we should take into account the scale dependence of the renormalized quantities a and 
g. The expression for the running coupling has the form 

(34) 

In the Lorentz gauge, the scale dependence of the gauge parameter a (in other words, the 
dependence of the normalization point) is governed by the normalization condition as follows: 
the propagator at = fi^ must have the form 



Gt.'iP) = ^ ((9.^ - ^ V ) . (35) 



From the equations for the characteristic curves and the expression (jH^ for the running 
coupling we obtain 

/,2\(13/22) 

3^1 + Cln(^) j 
where the constant C is determined from the normalization condition ()35|): 

C- , (37) 

p2\ 22 



(13 - 3ao) In 

Thus the scale dependence of the gauge-fixing parameter is given by 



cM^ ,, .. . Sjm^\ . (38) 



13-3ao(l-5)' \\n{pyA 

The dependence of the propagator on the normalization point and the momentum is obtained 
by substituting the expressions for the running gauge-fixing parameter ()38|) and running 
coupling (jH^ into formula After such substitution, we arrive at 

1 39 
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5.2 ihe case or maximal Abelian gauge 

As contrasted from the case of Lorentz gauge, the gluon propagator in the maximal Abehan 
gauge cannot be represented in the form (fT7|) : it can be parametrized by three functions as 
follows: Gf,Gf, and G^^). 



(40) 



^» = ?^f^--^VT^ + «^^f )• (41) 



These functions are not independent because there should be relations between them 
resulting from the Slavnov-Taylor identities. 

The renormalization group equation for the one-particle irreducible Green's function has 
the form 

r\ r\ r\ \ 

^'di^ + ^^^^d~g + ^"^^^9^ + ^^^^^^ + ^7 ^' = °' ^^^^ 

where G2 is one of the functions G°l^ , G*^^ , or "* and the renormalization group functions 
are defined by the equations 

m = ^, l{9) = ^ (43) 
2 9a 2^/3 



In the one-loop approximation, we obtain 



m = -hg\ (44) 



11 9' 

3 167r2' 



^o//^^) ^ „off^2 (17 -3a -6/?) 



Sa{g) 



12 167r2' 
(3a2-4a + 9) g^ 



IGtt 



2' 



The set of simultaneous equations defining the characteristic curves for the equation 
can be represented in a symmetric form, 

dX dg da dj3 dG2 



-1 I3{g) 5„ 50 -27G2 
Now we solve these equations for the Green's function of the "diagonal" gluon. Integra- 
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tion gives 



h 
h 
h 



2b\ + -, 
9^ 



arctan 



3a - 2 



11 



(3) 



Thus the general solution of the equation (jl^ has the form 



60 
40 
20 

-20 
-40 
-60 





; ' ' 1, = - — 

i ii = i 

I3-100 







10 20 30 40 50 60 70 80 

g 



(45) 



Figure 4: Projection of characteristic curves onto the {<y,g) plane. 



(46) 



where '^{x,y,z) is an arbitrary function. The solution of the equation (jl^ for is the 
same. A particular solution is obtained by a determination of the specific form of the function 
^(x, ?/,^) from the boundary condition. It is convenient to specify the boundary conditions 
on the plane = jj? {\ = 1) assuming that the perturbative expansion for G'^ is valid at 
small g. In the tree approximation, Grp' = 1 and, therefore. 



1 

yz' 



From this formula it follows that 



G 



(3) 



1 



(47) 



(48) 



3 167r2 /i2 

Taking the scale dependence of fjH^ into account yields 

(3) ^ In (^VA^) 
^ ln(pVA2)' 



(49) 
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The scale dependence of the parameter f3 can be determined from the equations for the 
characteristic curves and formula (jBll), the result is 



/3 



24 7r2 VA2 



(50) 



The parameter Cp is determined from the normalization condition /3 = /3o at = yU^, L-+ 



/3 = /3o 



ln(/iVA^: 
ln(pVA2) 



(51) 



The expression for the propagator for the "diagonal" gluon in the maximal Abelian gauge 
has a simpler form than the analogous expression in the Lorentz gauge. It can be accounted 
for by the remaining U{1) symmetry in the maximal Abelian gauge, which implies the 
relation Zg = Za^^'^ between the renormalization factors. Therefore, the functions f3{g) and 



j^^^g) are connected with each other by the formula P{g) = gj 



{3)_ 

Now we solve the equation for the transverse component of the Green's function in 
the case of "off-diagonal" gluon. First integrals h, I2, h of the system are the same as 
in the previous case. The fourth integral is given by 

= Gf (^i5/22gM^3^2 _ 4^ + 9)- i_ (52) 
Thus the solution of the equation ()42p . which is valid to both Cif,^ and G°[^ , takes the form 



Gf = (^ie- * (Sa^ - 4a + 9)^/^$ (Is, hji] 



(53) 



where $(x, z) is an arbitrary function, whose specific form is governed by the boundary 
conditions. For the sake of convenience, we set the boundary conditions on the plane = 
/i^(A = 1). It is assumed that, at small values of g, the perturbation theory expansion for 
the function G^-^ is valid. In the tree approximation (that is, to zeroth order in g), we obtain 



Gf = 1 and 



Therefore, 



3 / \/23 , 

— cos X in z 

23 I 22 



(54) 



where 



23 



p 



22 



g2/3/l67r2 



-15/44 



. 3a-2 . ^ 
cos 4 H 1=^ sm t, 



23 



1/2 



InC, C 



1 ^ ^^3' , 



(55) 



(56) 



The fi dependence of the parameters a and (3 is determined by the scale dependence of 
the running coupling (we restrict out attention to the one-loop approximation) and the 
values of the integration constants I2 and I3 in formulas ()45|) . A specification of these 
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values is equivalent to specification of the normalization conditions, which, in the case under 
consideration, is as follows: the propagator at = /i^ has the form 



Gab 



1 



In view of the relation (jBH), we obtain the equation 

ln(pVA2) 



p^ 



c 



ln(/iVA2) 



(57) 



(5^ 



giving the dependence of ^ on p^. Though it is a straightforward matter to derive the 
dependence of (3 on p^, the p^ dependence of a is derived in a more intricate way. To do this. 



we consider the factor 
Then we obtain 



. 3a-2 . ^ 
cos 4 H 1=^ sm 4 



23 



1/2 



and introduce the angle (p = arctan ■ 



3a -2 



cos 4 H 1=^ sm 4 



23 





cos(^ - ^) 




COS ip 



23 
(59) 



23 



From formula ()45|) it follows that p = + — — In (I^ is the integration constant 

independent of fi). Taking the scale dependence of the coupling constant into account, we 
arrive at 

V^(/^^)=/3-^ln(^ln(/iVA^)), (60) 

where the constant Is is determined from the normalization condition at = /x^: 

/ 11 



ip{p 



22 



In 



V247r2 



ln(pVA2^ 



Thus we obtain 



23 ^^/ ln(pVA^) \ 



(61) 



(62) 



22 Vln(/uVA2) ) 

and ^ — ^ = — v^o- From there formulas we derive the asymptotic behavior of the "off- 
diagonal" propagator: 

1/2 

I .'-I I \/ /j\ 1 t\\\\l}ll\\\ \ 

Qoff ^ ( PI] i^lnVMj 



^/A2) 
Vln(/iVA2) 



15 

44 



COS! ^In (MeIIA^ 

cos I 22 ill l^ln(MVA2) 



cosyjo 



(63) 



2 

The parameter (pQ in this formula is determined by the relation (po = arctan — ^_ , where 
ao is the gauge- fixing parameter at p"^ = jjp' (see (IHTjl ). 
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6 Discussion and conclusions 

The obtained expression ()63|) is singular at the values of momentum given by 



m{k) 



(64) 



13 



where k is an integer. An infinite number of singular points emerges botli wlien k — »• 
—oo {m{k) A) and wlien k —>■ oo {m{k) —>■ oo). In tlie neigliborliood of a singular 

point, the propagator behaves as C^^ ~ — — Neither such behavior nor the 

— m'^{k) 

existence of the singularities of a propagator in the Euclidean domain has a simple physical 
interpretation. Such points are probably artifacts of the normalization condition used here, 
one- loop approximation, and the maximal Abelian gauge in the continuum limit. It should 
be noted that a mathematically rigorous (say, as in ^T]) consideration of the question of 
whether this normalization condition is admissible can be the subject of a separate study. It 
should be emphasized that the maximal Abelian gauge in the continuum limit, that is, the 
condition 

id, + tgAl)iAl + tAl) = (65) 

provides a relation between the fields of the order An ^ 1 with the fields of the order An ^ -, 

9 

which may cause difficulties in computations in perturbation theory. The emergence of the 

terms — in the interaction Lagrangian is only one of manifestations of this relation. Such 

a 

terms emerge in the gauge-fixing part of the Lagrangian in a quite natural way. 

For this reason, the maximal Abelian gauge is essentially nonperturbative even in the 
domain of the small values of coupling constant. This conclusion is consistent with the 
renormalization-group analysis of the propagator behavior at the one-loop level. All the 
above imphes that the asymptotic behavior of the gluon propagator evaluated in this gauge 
on the lattice |3j, where nonperturbative contributions are naturally taken into account, may 
substantially diverge from the results of the renormalization-group improved computations 
in the framework of perturbation theory performed in the present study. 

The propagator of the " diagonal" gluon evaluated on a lattice jH] in the domain of large 
momenta can be adequately parametrized by formula PH)) with an appropriate choice of the 
parameters A,/i and Pq. 

The behavior of the "off- diagonal" propagator jS] diverges substantially form our results. 
The reason is as follows: 

1. There is no one by one correspondence between lattice MaG and perturbative MaG 
gauge; 

2. Renormalization-group improved computations do not take into account nonperturba- 
tive contributions, which are naturally taken into account in the lattice calculations. 

However, a comprehensive comparison of the renormalization-group improved perturba- 
tive behavior of the propagator with the behavior obtained numerically in a lattice gauge 
theory must include a consideration of the case when the gauge-fixing parameter a equals 
zero. In this case, the quantization procedure in the framework of the path-integral approach 
involves introducing an auxiliary field. A treatment of the gluon propagator at a = may 
form the subject of a subsequent study. 
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